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Abstract 

We study high energy cosmic rays from the decay of the gravitino dark mat- 
ter in the framework of supersymmetric model with i?-parity violation. Even 
though i?-parity is violated, the lifetime of the gravitino, which is assumed to 
be the lightest superparticle, can be longer than the present age of the uni- 
verse if i?-parity violating interactions are weak enough. We have performed a 
detailed calculation of the fluxes of gamma ray and positron from the decay of 
the gravitino dark matter. We also discuss the implication of such a scenario 
to present and future observations of high energy cosmic rays. In particular, 
we show that the excess of the gamma-ray flux observed by EGRET and the 
large positron fraction observed by HEAT can be simultaneously explained by 
the cosmic rays from the decay of the gravitino dark matter. 
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1 Introduction 



In particle cosmology, origin of dark matter of the universe is one of the most im- 
portant problems. Since there is no viable candidate for dark matter in the particle 
content of the standard model, new physics beyond the standard model is neces- 
sary to solve this problem. Supersymmetric model is a prominent candidate for the 
physics beyond the standard model; it not only introduces a viable candidate for 
dark matter, which is the lightest superparticle (LSP), but also solves other serious 
problems in particle physics, like naturalness problem of the electro-weak symmetry 
breaking. 

In order to realize LSP dark matter, conservation of i?-parity is usually assumed. 
In [1], however, it was pointed out that LSP dark matter scenario may be realized 
even with i?-parity violation (RPV) if the LSP is the gravitino; even though gravitino 
LSP becomes unstable with RPV, its lifetime may be longer than the present age of 
the universe because the decay rate of the gravitino is suppressed by the Planck mass 
as well as by (small) RPV parameter. Such a scenario has a great advantage for the 
thermal leptogenesis scenario [2j, as we will briefly discuss in the next section. Then, 
the primordial gravitino produced in the early universe can be a viable candidate for 
dark matter. 

Even though the lifetime of the gravitino is much longer than the present age 
of the universe, a fraction of the gravitinos have decayed until today. Such a decay 
becomes a source of high energy cosmic rays [U [3] ; the decay of the gravitino dark 
matter may produce high energy gamma ray and positron, which may be observed 
by present and future experiments. 

In this paper, we investigate gamma-ray and positron fluxes from the decay of 
gravitino dark matter in supersymmetric model with RPV. For this purpose, we first 
calculate the decay rate and branching ratios of the gravitino taking into account 
all the relevant operators. Fragmentation and hadronization of the decay products 
are studied by using PYTHIA package [4J. Then, we calculate fluxes of gamma ray 
and positron from the decay of the gravitino dark matter, carefully taking account 
of the propagation of the cosmic rays. We discuss the implications of our results 
to the present and future observations of the high energy cosmic rays. Importantly, 
excesses of the gamma-ray and positron fluxes over the backgrounds are reported 
by Energetic Gamma Ray Experiment Telescope (EGRET) [S] and High Energy 
Antimatter Telescope (HEAT) [6] experiments, respectively. We will show that these 
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excesses may be simultaneously explained by the scenario mentioned above. 

There are many past works to address HEAT and EGRET anomalies in terms of 
the annihilation of weakly-interacting-massive-particle (WIMP) dark matter [71 [H 
[9l [ini [11]. However, the predictions of the present scenario significantly differ from 
those of the annihilation scenarios. First, it is difficult to explain these anomalies 
simultaneously in the annihilation scenario; it is quite unlikely that dark matter 
annihilation is a main constituent of extragalactic gamma ray without exceeding 
the observed gamma-ray flux from the Galactic center pr2J. On the other hand, in 
the case of the decaying dark matter, the above constraint is relaxed, because the 
production rate of the gamma-ray is proportional not to the square of the number 
density of the dark matter but to the density. Then, as we will show, it is also 
possible to explain HEAT and EGRET anomalies simultaneously. It is also notable 
that, in the decaying dark matter scenario, the fluxes of the high energy cosmic rays 
are insentive to the boost factor, on which the fluxes in the annihilation scenarios 
strongly depend on. 

The organization of this paper is as follows. In Section [21 we first summarize the 
cosmological scenario that we consider. In Section [31 we discuss the decay processes 
of the gravitino. In Section [H formulae to calculate the cosmic-ray fluxes are given. 
The gamma-ray and positron fluxes from the decay of the gravitino dark matter are 
shown in Section [S] readers who are mainly interested in the results may directly go 
to this section. Section [SI is devoted to conclusions and discussion. 

2 Cosmological Scenario 

We first introduce the cosmological scenario that we consider. Although i?-parity 
conservation is usually assumed in conventional studies of supersymmetric models, 
RPV has a favourable aspect in cosmology. In supersymmetric models, it is often 
the case that the thermal leptogenesis scenario [2], which is one of the most promi- 
nent scenario to generate the present baryon asymmetry of the universe, is hardly 
realized since such a scenario requires relatively high reheating temperature after 
inflation, Tr > 10^ GeV [T3l [Ti]. With such a high reheating temperature, gravitino 
overproduction problem arises for wide range of the gravitino mass as far as i?-parity 
is conserved [12]. If the gravitino is unstable, gravitino produced after the reheating 
decays after the big-bang nulceosynthesis (BBN) starts and spoils the success of the 
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BBN 
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If the gravitino is stable, on the contrary, the primordial gravitino survives 
until today and contributes to the present energy density of the universe. Since 
the gravitino abundance increases as the gravitino mass becomes smaller, overclo- 
sure of the universe happens unless the gravitino mass is large enough [18]. With 
Tr > 10^ GeV, the above problems can be avoided only when (i) the gravitino is 
stable, and (ii) the gravitino mass is around 100 GeV. However, even in such a 
case, one has to worry about the decay of the lightest superparticle in the minimal 
supersjTiimetric standard model (MSSM) sector, which we call MSSM-LSP; with im- 
parity conservation, the MSSM-LSP, which is assumed to be the next-to-the-lightest 
superparticle (NLSP), decays only into gravitino and some standard-model parti- 
cle(s). When m3/2 ~ 100 GeV, the lifetime of the MSSM-LSP becomes longer than 
1 sec and relic MSSM-LSP decays after the BBN starts. When the MSSM-LSP is 
the neutralino or charged slepton, such decay processes spoil the success of the BBN, 
and hence it is difficult to realize the thermal leptogenesis scenario. 

If the i?-parity is violated, the MSSM-LSP may decay via RPV interaction and 
its lifetime may become shorter than 1 sec. Then, Tr ~ 10^ GeV is allowed. In such 
a case, gravitino is no longer stable and decays to standard-model particles. Even 
in such a case, however, the lifetime of the gravitino may be longer than the present 
age of the universe and hence the gravitino dark matter and thermal leptogenesis 
may be simultaneously realized [1]. 

Here, we consider the case where the gravitino is the LSP in the framework of 
the i?-parity violated supersymmetric models. We assume that the present mass 
density of the gravitino is equal to the observed dark matter density so that the 
gravitino can play the role of dark matter. There are several possibilities of the 
origin of such a primordial gravitino: scattering processes of thermal particles [TS] 
or the decay of scalar condensations [T9jJ. We will not discuss in detail about the 
production mechanism of the primordial gravitino because the following arguments 
hold irrespective of the origin of the gravitino. In addition, we consider the case that 
the lifetime of the gravitino is much longer than the present age of the universe; the 
upper bounds on the size of RPV couplings will be discussed in the following section. 



"^^However, ~ 10^ GeV may be also allowed when the gravitino mass is larger than 0(10 TeV). 
Such a scenario may be realized in the class of anomaly-mediation model 117j . 
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3 Framework and Decay Rates 



In this section, we introduce the supersymmetric model that we consider. Then, we 
summarize the decay rates of the gravitino and the NLSP, which are important for 
our study. 



3.1 Model 

In this article, we concentrate on the case where the i?-parity violating interactions 
originate from bi-linear terms of Higgs and lepton doublet. In the original basis, the 
i?-parity violating interactions are assumed to be bi-linear terms in superpotential 
and supersymmetry (SUSY) breaking terms [20]. Without loss of generality, we 
can always eliminate the bi-linear i?-parity violating terms from the superpotential 
by the redefinition of the Higgs and lepton-doublet multiplets. Then, the mixing 
terms between the Higgsino and lepton doublets are eliminated from the fermion 
mass matrix. In the following, we work in such a basis. Then, the relevant i?-parity 
violating terms are only in the soft-SUSY breaking terms, which are given by 

^RPV = BiLiHu + m\^^LiH*^ + h.c, (3.1) 

where Lj is left-handed slepton doublet in i-th generation, while and are up- 
and down-type Higgs boson doublets, respectively. In the following, we study the 
phenomenological consequences of the i?-parity violating terms given in Eq. (13.11) . 

With these i?-parity violating terms, the vacuum expectation values (VEVs) of 
left-handed sneutrino fields z>j are generated as 



Bi sin B + m? cos (3 
(^.) = ^ V, (3.2) 



where v ~ 174 GeV is the VEV of standard-model-like Higgs boson, tan/? = 
and TTiy. is the mass of z/j. The VEVs of the sneutrinos play impor- 
tant role in the following analysis. We parametrize the VEVs of the sneutrinos as 

^^ ^ ^, (3.3) 

V 

and consider the case that <^ 1. 

One important constraint on the size of the i?-parity violation is from the neutrino 
masses. Here, we assume that the neutrino masses are mainly from some other inter- 
action, like the seesaw mechanism [21j or Dirac-type Yukawa interaction. However, 
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the VEVs of sneutrinos also generate neutrino masses; assuming the Majorana-type 
masses for neutrinos, the ij component of the mass matrix receives the contribution 
of 

|c--o P 

[Am,] . = mlKiKj V (3.4) 



where rriz is the Z-boson mass. In addition, Zino Z, which is the superpartner of the 
Z-boson, is related to the mass eigenstates of the neutralinos Xa (with mass fn^o^) 
as 

^ = E^^xoX°. (3.5) 



(We also define the coefficients for photino and Higgsinos, c^y^o, cgo^o, and c/^o^o, 
by replacing Z — > 7, and H^-) Assuming that the neutralino masses are close 
to the electro-weak scale so that the SUSY can be the solution to the fine-tuning 
problem of the electro-weak symmetry breaking, the correction to the neutrino mass 
matrix is estimated as 

|Am„|, ~ 10- eV x (^) (^) . (3.6) 

It indicates that the i?-parity induced neutrino mass does not exceed experimental 
bound of observed neutrino mass when Kj < 10"'' is satisfied. 

As we have mentioned, one of the important motivations to consider RPV is 
to relax the BBN constraints due to the decay of the MSSM-LSP. In a case that 
MSSM-LSP is Bino-like neutralino B, it decays in two-body processes, B Zui, 



Wli, and hvi. The decay rates of each mode are given 



^B^z.. = -^9hin'ewK^m^[l-3^^+2^^], (3.7) 
^B^wi. = ^9hm'e^K^m^^l-3^ + 2^], (3. 

where gz = ^/ gf + gi (with gi and g2 being the gauge coupling constants of the 
f/(l)y and SU{2)l gauge groups, respectively), 6y/ is the Weinberg angle, is 



'^^Here, we consider the case that the hghtest Higgs boson h is almost standard-model like, so 
that the Higgs mixing angle is given by the [3 parameter. 
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Bino-like neutralino mass, and mx {X = Z, W, h) is the mass of gauge or Higgs 
boson. Hence the hfetime is estimated as 



0.01 sec X (^) - (^)-, ,3.10) 

where 

(3.11) 



In another case that MSSM-LSP is right-handed stau tr, it decays in the processes, 
fn TUi. In such a case, the decay rate is given by 

^fR = ^gzsin'^Owi^'^ (-^) rrif^, (3.12) 

where m^o and rrif^ masses of the lightest neutralino and stau, respectively, and the 
lifetime is estimated as 

~0.3secx( -) ( ^ ) ( ^) . 3.13 

^ VlO-"/ V300 GeV/ V200 GeV/ ^ ' 

Therefore, the lifetime of the MSSM-LSP, which is the NLSP in this case, becomes 
shorter than ~ 1 sec if typically Ki > 10^^^ is satisfied. Thus, combined with the 
upper bound for from the neutrino mass, we focus on the parameter region: 

10-'' <^^^< 10-'. (3.14) 

Before closing this subsection, we comment on the effects of tri-linear i?-parity 
violating terms induced by the redefinition of the Higgs and lepton-doublet multi- 
plets. With the redefinition of and Lj to eliminate the bi-linear i?-parity violating 
terms from the superpotential, tri-linear i?-parity violating terms are induced. They 
are irrelevant for our following studies, but are constrained, in particular, from the 



wash-out of the baryon asymmetry of the universe 
i?-parity violating terms in the superpotential as 
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Let us denote the tri-linear 



W^RPV — KjkLkLiEj + X'-jj^L^QiDj, (3.15) 

where Lj, and Qi are left-handed lepton, quark doublets, while E'^ and are right- 
handed lepton, down-quark singlets, respectively. (Here, "hat" is for superfield.) 



'^'^In the present setup, baryon number is conserved, so the constraints from the nucleon decays 
are irrelevant. 
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Figure 1: Diagrams of gravitino decay. 



Then, in order not to wash out the baryon asymmetry of the universe, the couphng 
constants in the above superpotential are constrained as [22] 

A,,fc,A^,<10-^. (3.16) 

For example, if we assume that the size of i?-parity violating terms are 0{Ki) relative 
to the corresponding i?-parity conserving ones (which are obtained by replacing 
with Li), and that the size of the SUSY breaking parameters are typically of the 
order of the electro-weak scale, then the above constraint is consistent with the one 
obtained from the neutrino mass. 

3.2 Gravitino decay 

In the case with RPV, gravitino LSP is no longer stable and decays to standard- 
model particles with a finite lifetime [23]. Here, we will take a closer look at the 
gravitino decay. 

In the present scenario, gravitino mainly decays in the two-body decay processes 
shown in Fig. ^ ^ •yi'i, Zui, Wli, and hui. (Here and hereafter, denotes the 
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gravitino.) Decay widths of each process are given bjcl 

r 



1 ^^^'^3/2 ^2^, 



7' 



2 n2 



1287r M^y 
12871 Ml^ 



(3.17) 



F{m3/2,mz) + 



3v 



2m 



'3/2 



>G{m3/2,mz) 



(3.18) 



647r 



3t; 



/T7i 



3/2 



4 K^^m? 



'3/2 



2 



3847r M^j 



m-, — m 



+ sm 



h 



(3.19) 
(3.20) 



where Mpi ~ 2.4 x 10^^ GeV is the reduced Planck mass, m3/2 is the gravitino mass, 



m 



2 ' 

3/2 



(3.21) 



and the functions F, G, and H are given by 



F{m3/2,mx 

H {7713/2, mx 

In addition, we define 

9^ : 



1 2 T 4 

1 mj^ 1 ^T^x 



1 m' 



X 



3 rng^g 3 m3^2 
^ 1 1 m\ 



o 6 ' 
.3 "^3/2 



2 m3^2 2 m3/2 



1 + 10 



m 



X 



4 



2 ' 4 
■"^3/2 "^3/2 



+ 



(3.22) 

(3.23) 
(3.24) 



m^o 

a=l 



a=l 
4 



1 NT^ '^W'+xJ'^W^-Xa 



+ h.c. 



(3.25) 
(3.26) 
(3.27) 



where c^± -± is the elements of unitary matrices which diagonalize the mass matrix 
of charginos M^±: = Tl=iCw±xiXt (i-e-, m^± = Eij Ci^jc.-± [A^^i]^^). 
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Figure 2: Lifetime of gravitino as a function of gravitino mass. Here, we take tan/3 = 
10, TJih = 115 GeV, nij^ = 1.5m3/2, rrii^ = 2m-^/2 under large Higgsino-mass limit, and 
assume GUT relation among gaugino masses. 

Lifetime of gravitino is determined by these two-body decay processes: 

3 

'^372 = ^3/2 = 2 ^ [r^^^^,,^ + T^^^zu, + T^^^wk + r^^hz/J , (3.28) 

i=l 

where the factor of 2 is for CP-conjugated final states. In Fig. [2], we plot as 
a function of ^3/2 in the large Higgsino-mass limit. Here, we take tan/5 = 10, 
rrih = 115 GeV, mc = 2m3/2, = 1.5m3/2, and grand unified theory (GUT) 
relation among the gaugino masses are assumed. One can see that r3/2 is much 
longer than the age of the universe (^ 4.3 x 10^'^ sec) for weak-scale gravitino mass 
when K < 10~^. Thus, in such a parameter region, most of gravitinos produced in 
the early universe survive until the present epoch. 

Even though the lifetime of the gravitino is long enough to realize the gravitino 
dark matter scenario, it may be possible to observe the decay of gravitino dark matter 
at present epoch. In particular, high energy photons and positrons are emitted in 
the decay processes as well as in the following cascade decay processes. To see this, 

"^^In. [5], the coupling of the gravitino to the supercurrcnt of the slepton multiplet, which gives 
rise to the terms proportional to the functions G and H, was neglected. Consequently, the decay 

rates of the gravitino into W^P and Zi' are underestimated, resulting in enhanced branching ratio 

for the process V'/j ~* l^^- 




1x13/2 (GeV) 
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Figure 3: Branching ratio for each decay mode. Lines with the indices 'W/," "Zi^," 
"/iz^" and "71/" show Br{^pf, W+1-) + 5r(?/;^ ^ W-1+), Br{ij^ Zu) + Br{^f, 
Zv), Br{ip^ hv) + Br{'ip^ — > hu), and Br{ilj^ — > 71/) + Br{ilj^ — > 71/), respectively. 
(Summation over the generation index is implicit.) Here, we take the MSSM parameters 
used in Fig. [2] 



in Fig. [31 we show branching ratio of each decay mode for tan /5 = 10 in the large 
Higgsino mass hmit. (Notice that the branching ratios are independent of Kj.) 

When m3/2 ^ 80 GeV, the decay mode ip^ 7Z/ dominates in total decay 
rate, because the decay processes with the emission of the weak or Higgs boson are 
kinematically blocked. On the contrary, once the gravitino becomes heavier than 
the weak bosons, the branching ratio for the process ip^j, 7Z/ is suppressed. This 
behavior can be understood from the fact that the 6'-parameters defined in Eqs. 
03.251) — fl3.27p . in particular, Ox^, is suppressed when the neutralino and chargino 
masses become much larger than the electro- weak scale. 

This fact has important implication in the study of the scenario using cosmic 
rays. When the gravitino is lighter than ~ 100 GeV or so, line spectrum of the 
gamma ray may be a striking signal. With larger gravitino mass, on the contrary, 
it becomes difficult to observe the gamma-ray line spectrum. Even in that case, 
however, significant amount of gamma rays with the energy of 0(1 — 100 GeV) are 
emitted in the cascade decays of W, and h bosons. Therefore, the continuous 
high energy gamma ray could be used as another characteristic signal. In addition to 
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the gamma ray, high energy positrons are also emitted, which is another interesting 
signal from the decay of gravitino dark matter. 



4 Cosmic-Ray Fluxes: Formulae 

As we have discussed in the previous section, energetic gamma and positron are 
produced if gravitino decays via i?-parity violating interactions. If gravitino is dark 
matter of the universe, such decay products can be a source of high energy cosmic 
rays. In order to discuss how well the scenario is tested by the use of cosmic ray, it 
is necessary to formulate the calculation of the cosmic-ray fluxes from the decay of 
gravitino. In this section, we show how we calculate the gamma-ray and positron 
fluxes from the gravitino decay. 



4.1 Gamma ray from the gravitino decay 



The total flux of the gamma ray from the decay of dark matter (i.e., gravitino) is 
calculated by the sum of two contributions: 



dE 



DM 



dE 



+ 



dE 



(4.1) 



halo 



where the first and second terms in the right-hand side are fluxes of gamma ray from 
cosmological distance and that from the Milky Way halo, respectively. We discuss 
these contributions separately. 

The flux of the gamma ray from cosmological distance is estimated as 

E^ r ,^dN^{E') 



dE 



dE'G^i^E , E' 



dE' 



(4.2) 



Here, the propagation function of gamma ray is given by 



G,iE,E') 



Cpc^3/2 1 f E\ 



A-kH^VLmE \E'J 



3/2 



'i + njiiMiE/E') 



(4.3) 



where c is the speed of light, Hq is present Hubble expansion rate, pc is critical 
density, and ^3/2 ~ O.lUSh-^, Qm ^ 0.1369/1"^ ~ 0.721 (with h ~ 0.701) 
are density parameters of gravitino dark matter, total matter, and dark energy, 
respectively [21] • In addition, dN^/dE is the energy spectrum of gamma ray from 
the decay of single gravitino, which can be given by the sum of contributions from 
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relevant decay modes: 



'dE 






'dN^' 


r 





dN^ 



+ r 



+ r. 



dN^ 

He 



'dN^' 






hui j 



(4.4) 



Here, [diV^/d-E'] are energy distributions for each decay modes. Notice that the 
dN^/dE is determined once the SUSY parameters are fixed, irrespective of the cos- 
mological scenario. We calculate dN^/dE by using PYTHIA package [1]. 

The fiux of the gamma ray from the Milky Way Galaxy halo is obtained as [25] 

E^ 



E' 



) dJ -y 

liE 



1 dN^ 
47r~dE 



l.o.s. 



Pz/2il)dl ) 



(4.5) 



/ dir 



halo "^3/2^3/2' 

where P3/2 is the energy density of the gravitino in the Milky Way halo. In Eq. (14.51) . 
the integration should be understood to extend over the line of sight (l.o.s.). Thus, 
the integration has an angular dependence on the direction of observation. Here, 
[E^dJ-y/dE]iiaio is given by averaging over the direction, which is denoted as (■ ■ • )dir- 
In the EGRET observation, the signal from the Galactic disc is excluded in order 
to avoid the noise. In order to compare our results with the EGRET results, we 
also exclude the region within ±10° around the Galactic disk in averaging over the 
direction. 

In order to perform the line-of-sight integration, the profile of P3/2, namely dark 
matter mass density profile Phaio? should be given. For our numerical analysis, we 
adopt Navarro- Frenk- White (NEW) density profile [26]: 

Ph 



(4.6) 



r/rc(l + r/r^)^ ' 

where r is the distance from the Galactic center, ph — 0.33 GeV cm~'^, and Vc — 20 
kpc. We have checked that the dependence on the dark matter profile is negligible 
because, in the calculation of the gamma-ray fiux, we exclude the region around the 
Galactic disc as we have mentioned. 



4.2 Positron from the gravitino decay 

Next, we discuss cosmic-ray positron from the gravitino decay. If we consider ener- 
getic positron propagating in galaxy, its trajectory is twisted because of magnetic 
field. With the expected strength of the magnetic field, scale of gyro-radius of the 
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trajectory is much smaller than the size of the galaxy. Furthermore, the magnetic 
field in the galaxy is entangled. Because of these, propagation of the positron in the 
galaxy is expected to be well approximated as a random walk. 

We use a diffusion model for the propagation of positron, in which random walk 
is described by the following diffusion equation [SI [TT] : 

^^^^ = K{E)V'U{E,x) + ^ [b{E)UiE,x)] + QiE,x), (4.7) 

where fe+{E,x) is the number density of positrons per unit energy (with E being 
the energy of positron), K{E) is the diffusion coefficient, b{E) is the energy loss 
rate, and Q{E,x) is the positron source term. As we mentioned above, diffusion of 
injected positron is caused by the entangled magnetic field in the galaxy. On the 
other hand, the energy loss of the positron is via Thomson and inverse Compton 
scatterings with Cosmic Microwave Background and infrared gamma ray from stars 
or Synchrotron radiation under the magnetic field. The functions K{E) and b{E) 
can be determined so that the cosmic-ray Boron to Carbon ratio and sub-Fe to Fe 
ratio are reproduced. In our analysis, we use those given in [8]: 



K{E) = 3.3 X 10^^ 



1.39 



E 



1 GeV 



cm sec , (4. 



Since the magnitude of the energy loss rate indicates that positron loses its energy 
in the fiight of less than a few kpc, the positron fiux from outside of our Milky 
Way Galaxy halo is negligible. Thus, in the following discussion, we focus on the 
contribution of the positron fiux from the Milky Way Galaxy. In addition, the 
positron source term is given by the use of the positron injection rate and dark 



matter distribution in the Milky Way Galaxy halo aa^^ 



QiE,x) = (^^^^, (4.10) 

m3/2 Ty2 dE 

where dNe+ / dE is energy distribution of positron from the decay of single gravitino. 



*^The uncertainty on the positron flux from the effect of inhomogeneity in the local dark mat- 
ter distribution is negligible [27^, which is quite contrast to the traditional case of dark matter 
annihilation. 
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The explicit expression is given as 



dNe+ 
~dE 




J w-v 



(4.11^ 



where [dNf,+ / dE\... is energy distribution for each decay mode. We calculate the 
energy distributions by the use of PYTHIA package. For density distribution, we 
adopt the same profile as the previous section, namely NFW profile defined in Eq. 

We solve the diffusion equation in finite diffusion zone with boundary condition. 
Since the observed cosmic-ray positrons are considered to be in equilibrium, we 
impose stability condition df{E, x)/dt = 0, and also free escape condition f{E, x) = 
at the boundary. The diffusion zone is usually assumed as cylinder characterized 
with half-height L and radius R. Since positrons lose their energy after the flight of 
a few kpc or less, the positron flux does not strongly depend on the choice of the 
diffusion zone. In our analysis, we take L = 4 kpc and R = 20 kpc. 

Positron flux from the decay of gravitino dark matter is given by 

dJp^ 



DM 



dE 



-^fiE,Re), 



(4.12) 



where Rq is the location of the solar system. This flux does not correspond exactly 
to the one observed on the top of the atmosphere. The flux is modified due to 
interaction with solar wind and magneto-sphere. However, the modulation effect is 
not important when the energy of a positron is above 10 GeV. Furthermore, the 
effect is highly suppressed in the positron fraction, which is defined by the ratio of 
the positron flux to the sum of positron and electron fluxes, i.e., $e+/('^'e+ + '^'e-)- 



4.3 Backgrounds 

In order to discuss high energy gamma ray and positron as signals from the decay 
of dark-matter gravitino, it is important to understand those cosmic rays from other 
sources. (We call them backgrounds.) 

Cosmic gamma rays have various origins. As we have mentioned, the cosmic 
gamma rays can be divided into two parts by their origins; Galactic and extragalactic 
parts. Some part of the Galactic origins, such as scattering processes (i.e., inverse 
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Compton scattering and bremsstrahlung) and pion decay, are known as probable 
sources of cosmic gamma rays. However, other Galactic origins, as well as the extra- 
galactic ones, have not been well understood. In the study of the gamma ray from 
the gravitino decay, it is important to understand the behavior of the unidentified 
cosmic gamma ray (UCGR) from various originsEj However, theoretical calcula- 
tion of such components of the cosmic gamma ray is difficult, so we adopt more 
phenomenological approach to extract the UCGR not originating from the gravitino 
decay. 

Currently, the cosmic gamma ray flux has been measured by EGRET, and various 
analysis have been performed to extract the UCGR from the EGRET data. The first 
intensive work was done in [5], in which it is concluded that the UCGR follows a 
power law as E^dJ^/dE = 1.37 x 10-'^(E/1 GeV)-^-\cm^ str sec)~iGeV. Generally, 
however, there is a difficulty in removing the contribution from the known scattering 
or decay processes in the Milky Way Galaxy since the analysis depends on the 
Galactic model. Recently, with an improved analysis in the estimation of the Galactic 
contribution [33l it has been pointed out that the UCGR spectrum follows a 
power-law in the energy range E < 1 GeV. However, for > 1 GeV, a deviation 
from the power-law behavior is reported. As we will see, the effect of the gravitino 
decay on the cosmic gamma ray becomes important for the energy range of E > 
0.1 — 1 GeV. Thus, in our analysis, we assume that the UCGR in the lower energy 
range is only from astrophysical origins (although many of them have not yet been 
well understood). In addition, we also assume that the spectrum of UCGR from 
astrophysical origins follows a power law, and hence its behavior can be extracted 
from the data in the sub-GeV region. Since the gamma ray from the gravitino 
becomes important above the energy of 0.1 — 1 GeV, we use the observed data in the 
range of 0.05 GeV < < 0.15 GeV to determine the background flux. Assuming 



the power-law behavior, we obtain the best-flt UCGR flux as 

/ \ -0-449 

~ 5.18 X 10"^ (cm^ sec str)"^ GeV x \ ■ (4-13) 



2 dJ'y 

^ lE 



"^^The UCGR may have various origins. Examples in astrophysics are contributions from galaxy 
clusters [5S] , energetic particles in the shock waves associated with large-scale cosmological structure 
formation [29j . distant gamma-ray burst events, baryon-antibaryon annihilation |30j . In addition, 
if we consider physics beyond the standard model, spectrum of UCGR may be affected by, for 
example, the evaporation of primordial black holes [31], the annihilation of WIMPs [7l [8l [9l ITOl fTT] . 
or extragalactic IR and optical photon spectra In our analysis, we only consider the decay of 
the gravitino as a particle-physics source of the UCGR, and do not consider other possibilities. 
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We use this spectrum as the background in the following analysis. Then, the total 
gamma-ray spectrum is given by 









+ 


dJ^ 


_dE_ 


tot 


_dE_ 


DM 


_dE_ 



(4.14) 



BG 



Next, let us consider the background positron (and electron). Cosmic rays mainly 
consist of nuclei, electrons, and positrons. Nuclei are the dominant component of 
the cosmic rays and pouring to the earth after it has drifted by interaction with 
interstellar matters in our Galaxy. As a consequence, secondary cosmic-ray electrons 
and positrons are produced. On the theoretical side, many simulations of cosmic-ray 
electron and positron have been done by the use of cosmic-ray propagation model 
|33] . In our study, we adopt the following cosmic-ray electron and positron from 
astrophysical processes [8]: 



J prim 



l + llKO-9 



GeV 



+ 3.2^2.15 



(GeV cm^ sec str) ^, 



(4.15) 



GeV 



1 + llOE^^V + 600^2 ev + 580E^g!y 
1 + 650^2-3^ + 1500^4-2^ 



(GeV cm2 sec str)-\ (4.16) 
(GeV cm2 sec str)~\ (4.17) 



where -EgbV is the energy of electron or positron in units of GeV. We note here that 
concerning the backgrounds, the secondary electron accounts for about 10 % of the 
total electron flux while the positron flux is dominated by the secondary one. 

With these backgrounds, the total fluxes of the electron and positron are obtained 

as 



['^e+]DM + ['^'e+]sec. 
['^e-]DM+['^e-]pHm+['^. 



(4.18) 

(4.19) 



Importantly, the above background in cosmic-ray positron flux well agrees with 
HEAT observation in the energy range E < 10 GeV. However, for E > 10 GeV, an 
excess of the positron flux is seen in the HEAT data [B]. In the following section, we 
show that such an excess may be due to the decay of gravitino0 



#'^For other possibilities, see [71 [8l [9l ITT] . 
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5 Numerical Results 



5.1 Cosmic rays from the gravitino decay 

Now we are at the position to present our numerical results. In this subsection, we 
show fluxes of cosmic-ray gamma and positron from the gravitino decay to discuss 
their behaviors. 

First, we discuss the gamma-ray flux. The flux for m^i2 = 150 GeV is shown on 
the left in Fig. HI (Here, the lifetime of t^i2 = 10^^ sec is used.) We found that the 
dependence of the gamma-ray flux on the flavor of the primary lepton is negligible. 
(Here and hereafter, for the calculation of the gamma-ray flux, we assume that the 
gravitino mainly decays into third-generation leptons and gauge or Higgs bosons.) 
In calculating the gamma-ray flux, we adopt the energy resolution of 15 %, following 
EGRET [5] . In the figures, we also show the contributions of individual decay modes. 
(In this case, the contribution of the decay mode into hv is very small.) 

As one can see, continuous spectrum is obtained from the decay modes into weak 
boson (or Higgs boson) and lepton, while a relatively steep peak is also obtained at 
E = |w^3/2 due to the monochromatic gamma emission via 'ju. With larger 

gravitino mass, Br{ip^ — > 7//) is suppressed, and hence the peak a.t E = \m^/2 
becomes less significant. Thus, when the gravitino mass becomes much larger than 
the masses of weak bosons, it will become difficult to find such a line spectrum. 

The positron fraction is shown on the right in Fig. H] for the cases where the 
gravitino dominantly decays into first, second, and third generation leptons. Since 
the energy of positron becomes smaller during the propagation in the Galaxy, the 
positron spectrum has the upper end point at ~ |w^3/2- In particular, when the 
gravitino can directly decay into positron, the end point becomes a steep edge; such 
an edge can be a striking signal of the decaying dark matter if observed. Even in 
other cases, the positron spectrum has a peak just below the upper end point. This 
may also provide an interesting signal in the observed positron flux. 

As one can see, the fluxes of the gamma-ray and positon in the present scenario 
are most enhanced for the energy of 0(1 — 100 GeV). In addition, the fluxes are 
significantly suppressed for the energy smaller than ~ 1 GeV. Thus, the present 
scenario is not constrained from the observations in such a low-energy region, like 
the observations of the fluxes of sub-GeV gamma rays. 

So far, we have shown results only for the case of r3/2 = 10^^ sec. Fluxes with 
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0.1 1 10 100 1000 0.1 1 10 100 1000 

E (GeV) E (GeV) 



Figure 4: Left: Gamma-ray flux from gravitino decay [£'^(iJ-y/(i£']DM("total"). Lines with 
"71/," "Zz/," and 'Wr' are contributions from each decay mode. Right: Positron fraction 
[$e+]DM/([^e+]tot + [^^'e-Jtot)- Lines with "e," and "r" show the results for the case 
that the gravitino mainly decays to first, second, and third generation leptons, respectively. 
For both of the figures, we take = 150 GeV, r3/2 = 1.0 x 10^^ sec, and the MSSM 

parameters used in Fig. [21 

other values of the lifetime can be easily obtained from the above results since the 
cosmic-ray fluxes from the gravitino decay is inversely proportional to One 
important point is that, when r3/2 = 0(10^^ sec), both the gamma-ray and positron 
fluxes from the gravitino decay become comparable to the background fluxes dis- 
cussed in the previous section for ~ 1 — 100 GeV. Thus, with such a lifetime, we 
may be able to see the signal of the decay of gravitino dark matter in the spectrum 
of cosmic rays. The following subsections are devoted to the discussion of such an 
issue. 

5.2 Implications to present observations 

As discussed in the previous section, some anomalies are indicated both in the 
gamma-ray spectrum observed by EGRET and the positron fraction observed by 
HEAT. In this subsection, we show that these anomalies may be simultaneously 
explained by a single scenario, the gravitino dark matter scenario with RPV. 

In Fig. [5|, we show the total gamma-ray flux and the positron fraction for m^j2 = 
150 GeV and r3/2 = 2.2 x 10^^ sec. (We have used the best fit value of r3/2 for the 
EGRET data.) Here, we consider simple cases where the gravitino decays only into 
one of the three lepton flavors; only one of the (z/j) (i = 1 — 3) is non-vanishing 
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while the others are set to be zero. (For gamma ray, the total flux is independent 
of generation indices as we have mentioned in the previous subsection.) With the 
lifetime adopted, the gamma-ray flux and the positron fraction both signiflcantly 
deviate from the background. In gamma-ray flux (Fig. [5], left), one can see that the 
continuous spectrum originating from the processes ipfi ~^ and Zu gives a good 
agreement with EGRET data for E ^ 1 — lOGeV. In the positron fraction (Fig. 0, 
right), the results indicate that clear signal can be seen in the energy region E >10 
GeV over the background for all the three cases. In the same flgure, we also show 
the observational data of EGRET or HEAT. As one can see, agreements between the 
theoretical predictions and observations are improved both for EGRET and HEAT. 

Results for other values of the gravitino mass (and lifetime) are shown in Figs. 
[6] and [71 We can see that the suggested anomalies in the gamma-ray and positron 
fluxes may be explained in a wide range of the gravitino mass. 

In order to see the preferred parameter region in the light of EGRET result, we 
calculate the variable as a function of 771^/2 and 



where Xth,i is the theoretically calculated flux in i-th bin, which is calculated with 
Eq. (14.141) . Xobs,i is the observed flux, and (Tobs.i is the error of Xohs,i- In addition, 
is number of bins; A^ = 10 for EGRET. In Fig. [HI we show the region with < 18.3 
on the m3/2 vs. plane, which is 95 % C.L. allowed region. As one can see, the 
present scenario could well explain the EGRET anomaly in a wide parameter region, 
10^^ sec < Ts/2 ^ 10^'' sec and 7713/2 > 90 GeV. From Fig. [21 it can be seen that 
10~^° < Hi < 10^® is favored. In Fig. [HI we also show the parameter region which is 
consistent with the HEAT data (A^ = 9) at 95 % C.L. (i.e., < 16-9). 

As one can see, the present scenario can simultaneously explain the observed 
gamma and positron fluxes. 

5.3 Future prospects 

In the previous subsection, we have shown that the gamma-ray and positron fluxes 
from decaying gravitinos can successfully explain the results of the past observations. 
As we have seen, however, the energy ranges of the past observations are limited up 
to 0(10 GeV) although the signal from the gravitino decay may signiflcantly affect 
the cosmic-ray spectra up to the energy of 0(100 GeV). In addition, it is also true 




N 



(5.1) 
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Figure 5: Gamma-ray flux (left figure) and positron fraction (right figure). Here, we take 
m3/2 = 150 GeV, = 2.2 x lO^^ sec, and MSSM parameters as Fig. [H 
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Figure 6: Same as Fig. O except for 772,3/2 = 300 GeV and T3/2 = 1-5 x 10^^ sec. 



that the uncertainties of the cosmic-ray spectra observed by the past observations 
are relatively large at the energy range of ~ 0(10 GeV). Thus, it is desirable to 
test the scenario of gravitino dark matter with RPV with better observations. 

Fortunately, in the near future, new observations of cosmic rays. Gamma-ray 
Large Area Space Telescope (GLAST) and Payload for Antimatter Matter Explo- 
ration and Light nuclei Astrophysics (PAMELA), are expected to provide results 
of new measurements of the cosmic-ray fluxes. These experiments are designed to 
detect cosmic rays with energy up to a few hundreds GeV. Thus, they will give us 
better test of the scenario. Since GLAST (PAMELA) has better energy range and 
resolution than EGRET (HEAT) in the measurement of gamma-ray (positron) flux. 
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Figure 7: Same as Fig. El except for 777,3/2 = 500 GeV and r3/2 = 1-2 x 10 ° sec 



they should confirm the anomahes if they really exist. 

Even if the fluxes of the cosmic rays are smaller than the best-fit value of those 
observed by EGRET and HEAT, we still have a chance to see signals from the decay 
of dark-matter gravitino. To see expected constraints on the parameter space, we 
calculate the expectation value of the variable defined as 



N 



i=l 



BG,j 



(5.2) 



where is the number of bins. We use 50 bins to estimate the expected sensitiv- 
ities of GLAST and PAMELA. With the use of GLAST and PAMELA instrument 
parameters [53, EHl [ID], we define z-th bin as < E^j^+^i < where 



E. 



(min) 



7,1 



E 



(min) 



0.02 GeV X 



0.05 GeV X 



300 GeV \ 



GeV J 



0.02 



270 GeV \ 
0.05 GeV/ 



(5.3) 
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and E 



(max) 



E 



(min) 



^■y/e+ i ~ -^7/e+i+i- addition, Ath.i is the number of theoretically calculated 
events in i-th bin, which is the expected number of events in GLAST or PAMELA 
experiment. Notice that Ath,i is calculated by using Eq. fl4.14p or Eq. (14.181) . Fur- 
thermore, N^G,i is the number of background events. Here, we assume that the 
number of background events will be well understood in the future observations by 
using the date in low-energy range; in our following study, we use Eqs. (14.131) and 
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Figure 8: 95% C.L. allowed regions on vs. T3/2 plane; the shaded regions are the 

allowed regions. The region with "7" is for the EGRET data, while the region with "e" 
("^," "r") is for HEAT data for the case where the gravitino decays into first- (second-, 
third-) generation lepton. 
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(14.171) . Then, we obtain 

AE^^iTS, (5.5) 

tot 

AE^^iTS, (5.6) 

BG 

for gamma-ray flux, and 

iVth,. = [^e+{Ee+,^)],,,AE,+ ,TS, (5.7) 
iVBG,^ = [^eAEe+,^)iec^Ee+,^TS, (5.8) 

for cosmic-ray positron flux. Here, T and 5* are exposure time and acceptance, 
respectively, and ctbg,* is the error of the number of backgrounds, which we take 
c"BG,i = a/ Nbg,! assuming that the error is dominated by statistics. In our analysis, 
we take TS = 10^° and 10® cm^ sec str for GLAST and PAMELA, respectively. In 
addition, AE^j^+^i = E^^^\ — E^^^. is the width of the i-th bin. 

Expected constraints on the ^3/2 vs. plane is shown in Fig. [91 For positron 
flux, we consider the case that the primary lepton emitted by the gravitino decay 
is r or u^-. In Fig. [9l we show contours of (x^) = 67.5, corresponding to 95 % 
level of detectability. From the figure, one can see that GLAST and PAMELA 
have sensitivities for the case with the lifetime of 0(10^'' sec), and hence the signal 
from decaying gravitino can be observed in both observations in wide range of the 
parameter space. We have seen that, in order to explain the EGRET and HEAT 
anomalies, is required to be ~ 10^^ sec, which is order of magnitude shorter 
then the reach of GLAST and PAMELA experiments. Thus, even if EGRET or 
HEAT anomaly somehow dissapears, GLAST or PAMELA still has a chance to find 
some high-energy-cosmic-ray signals from the gravitino dark matter scenario with 
i?-partity violation]^ 

'^^Recently, the positron flux measured by PAMELA has been reported [37]. Importantly, the 
PAMELA result is still preliminary, and the positron flux is shown only up to the energy of 10 GeV. 
The new result is inconsistent with the previously calculated background flux, so the understanding 
of the background will become very important if the present result will be confirmed with more 
data. Thus, even though the PAMELA result is also inconsistent with the HEAT result, we think 
that it is still premature to discuss its implication. 
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Figure 9: Contour plot of (x^) on 771,3/2 vs. r3/2 plane for 95% probability to detect 
gamma-ray/positron signal from gravitino decay in GLAST/PAMELA. 



6 Conclusions and Discussion 

In this paper, we have discussed the high energy cosmic rays from the decay of 
gravitino dark matter in i?-parity violated supersymmetric model. Here, we have 
considered the case that the gravitino is the LSP, and that RPV originates from 
the Ljif„-type operators. In such a model, the gravitino dominantly decays into a 
gauge boson {W^, Z, or 7) and a lepton with extremely long lifetime. In particular, 
if RPV interactions are weak enough, the lifetime of the gravitino becomes much 
longer than the present age of the universe. Consequently, if the right amount of 
the gravitino is produced in the early universe, gravitino can be dark matter even 
if the -R-parity is not conserved. In addition, in such a scenario, the MSSM-LSP, 
which is assumed to be the NLSP, can decay with a lifetime shorter than ~ 1 sec via 
RPV operators. Then, the serious effects on the light element abundances by the 
decay of the NLSP, which gives one of the most stringent constraint on the scenario 
of gravitino dark matter, can be avoided. 

We have studied the gamma-ray and positron fluxes from the decay of the grav- 
itino dark matter. In our analysis, we have calculated the decay rate and branching 
ratios of the gravitino, taking account of all the relevant operators. Then, the energy 
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spectrum of the primary decay products (i.e., gauge bosons, leptons, and partons) 
are calculated. Decay and hadronization processes of those primary particles have 
been treated by using PYTHIA package to accurately calculate the primary fluxes of 
gamma ray and positron by the decay of gravitino. Then, by solving the propagation 
equations, we have obtained the fluxes of gamma ray and positron in the cosmic ray. 

One of our important results is that the anomalies observed by the EGRET and 
HEAT experiments can be simultaneously explained in this scenario if the lifetime of 
the gravitino is 0(10^^ sec). This conclusion holds for a wide range of the gravitino 
mass; as far as m^j2 is larger than ~ 80 GeV so that the gravitino can decay into 
weak boson(s), such a scenario works. In addition, with PAMELA and GLAST, 
more accurate test of the scenario will become possible; they will cover the parameter 
region of r3/2 ^ iQ^'i'-^s 

We comment here that, in the present scenario, the MSSM-LSP decays via RPV 
interactions, which may affect the LHC phenomenology. If the EGRET and HEAT 
anomalies are due to the decay of the gravitino dark matter, the lifetime of the 
MSSM-LSP is estimated to be 0(10^^ sec). This fact may have an impact on the 
LHC experiment. The typical decay length of the MSSM-LSP is 0(10'^ m), which is 
much longer than the size of the detector. Thus, most of the MSSM-LSPs produced 
at the LHC experiment escape from the detector before the decay. If the lightest neu- 
tralino is the LSP, the collider signatures are the same as the conventional signatures 
with the neutralino LSP. However, in the present scenario, the MSSM-LSP may be 
charged (or even colored). In such a case, we will observe a heavy charged particle 
as high pt track at the LHC. In addition, 0.1 — 1 % of the produced MSSM-LSP 
may decay in the detector, resulting in isolated vertices from the interaction point. 
These are very exotic signals which are not expected in the conventional model of 
supersymmetry. 

So far, we have considered the case where the gravitino is dark matter. However, 
with other candidates for dark matter, it may be possible to simultaneously explain 
the excesses of the gamma-ray and positron fluxes observed by EGRET and HEAT 
experiments. In particular, if dark matter decays mainly into the weak bosons with 
the lifetime of ~ 10^^ sec, this can be the case. One of the examples may be the 
lightest neutralino. With RPV, the lightest neutralino decays even if it is the LSP; 
with the RPV operator given in Eq. (13.11) . for example, the lifetime of the neutralino 
LSP becomes ~ 10^^ sec when k ~ 10~^^, assuming that the lightest neutralino is 
Bino-like and that the mass of the lightest neutralino is ~ 100 GeV. 
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Note added: While finalizing this paper, we found [38] which also studies high energy 
cosmic rays from the decay of the gravitino dark matter. In this paper, antiproton 
flux is also calculated, and it is discussed that the antiproton flux may become 
too large if we explain the EGRET and HEAT anomalies in the present scenario. 
However, it is also mentioned that the predicted antiproton flux suffers from large 
uncertainties. Thus, according to [M], taking into account the uncertainties, the 
antiproton flux can be consistent with the observations even if we explain the EGRET 
and HEAT anomalies. In addition, in the study of the gravitino decay in [38j, effects 
of the coupling of the gravitino to the supercurrent of the slepton multiplet have not 
been taken into account. 
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